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The interaction of a single manganese impurity with silicon is analyzed in a combined experimental 
and theoretical study of the electronic, magnetic, and structural properties of manganese-doped 
silicon clusters. The structural transition from exohedral to endohedral doping coincides with a 
quenching of high-spin states. For all geometric structures investigated, we find a similar dependence 
of the magnetic moment on the manganese coordination number and nearest neighbor distance. 
This observation can be generalized to manganese point defects in bulk silicon, whose magnetic 
moments fall within the observed magnetic-to-nonmagnetic transition, and which therefore react 
very sensitively to changes in the local geometry. The results indicate that high spin states in 
manganese-doped silicon could be stabilized by an appropriate lattice expansion. 

PACS numbers: 36.40.Cg, 75.50.Pp, 73.22.-f, 61.46.Bc 



The interaction of a deliberately introduced impurity 
with a semiconductor material is one of the most fun- 
damental problems of semiconductor physics. For mag- 
netic impurities, an important question is the survival or 
quenching of the magnetic moment. From the pioneer- 
ing studies of Ludwig and Woodbury ! 1J up to the present 
day, a wealth of experimental and theoretical studies have 
therefore been devoted to magnetic properties of transi- 
tion metal doped semiconductors [2~HB], semiconducting 
nanocrystals [7HH], and clusters [TUHIl]- A possible cor- 
relation between local magnetic moment and coordina- 
tion number of the impurity has been noticed in theo- 
retical work [BJ, but is challenging to investigate in bulk 
samples because of inhomogeneities, coalescence, or im- 
purity band formation. These difficulties can be over- 
come by utilizing size-selected, singly-doped clusters as 
model systems where a transition metal atom occupies 
a well-defined position in the silicon host, without any 
interaction between impurities. Here, we study MnSi^ 
by x-ray absorption and x-ray magnetic circular dichro- 
ism (XMCD) spectroscopy of size-selected free clusters 
[T5HT9] as a local and element-specific probe of electronic 
structure and magnetic properties. These experimental 
techniques are combined with non-empirical density func- 
tional theory (DFT) calculations. We find a clear depen- 
dence of the magnetic moment on the manganese coordi- 
nation and nearest-neighbor distance. This result can be 
generalized to manganese point defects in bulk silicon. 
Details of the experimental setup are given elsewhere 
[TBI I2U] . Very briefly, a continuous beam of MnSi^ clus- 



ters is produced in a magnetron gas aggregation source 
and transmitted through a combined radio-frequency 
hexapole ion guide and collision cell into a quadrupole 
mass filter. After mass selection, the clusters are ac- 
cumulated in a cryogenic linear Paul trap and thermal- 
ized to 10 — 20 K by collisions with helium buffer gas 
at p m 10~ 3 nrbar. To study the local electronic and 
magnetic properties of MnSi^ by x-ray absorption and 
XMCD spectroscopy, a tunable monochromatic x-ray 
beam delivered by an undulator beamline at the syn- 
chrotron radiation facility BESSY II is coupled on-axis 
into the ion trap for resonant excitation at the manganese 
L 2 ,3-edge. This creates Mn + and SiJ photoions, which 
are detected by a reflectron time-of-flight mass spectrom- 
eter. The incident photon energy is scanned from 618 - 
686 eV to record photoion yield spectra that are a mea- 
sure of the x-ray absorption cross section. For XMCD 
spectroscopy, which requires alignment of the total mag- 
netic moment of free MnSi^, the liquid- helium cooled 
ion trap is placed inside the homogeneous magnetic field 
(B = 5 T) of a superconducting solenoid, and ion yield 
spectra are recorded for parallel and antiparallel align- 
ment of photon helicity and magnetic field |21j . 
In addition to magnetic and electronic properties, struc- 
tural properties of MnSi^ are investigated. Similar to the 
reactivity and adsorption studies of doped silicon clus- 

I, the 



ters by Ohara et al. [22] and Janssens et al. 



exohcdral-to-endohedral transition of manganese-doped 
silicon cluster cations is monitored via the depletion 
of MnSi^ in the cluster beam when introducing p w 
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FIG. 1. (color online) Manganese 2p x-ray absorption (left) 
and XMCD (center) spectra of MnSijt clusters (n = 7 — 14), 
indicating quenched magnetic moments for n > 11; corre- 
sponding ground state structures of MnSi^ and the depletion 
of singly doped clusters in the presence of O2 as a measure of 
the exohedral-to-endohedral transition (right). 



10~ 3 mbar partial pressure of oxygen reactant gas into 
the hexapole collision cell. As can be seen in Fig. [l] 
the depletion of singly doped MnSi^ is 89 — 94 % for 
n = 7 - 10 but drops to - 15 % for n > 11. This is 
due to the large difference of manganese and silicon reac- 
tivity towards oxygen that makes this depletion study a 
highly sensitive measure of the exohedral-to-endohedral 
transition, which takes place from MnSif to MnSi^. 
This structural transition coincides with a marked change 
in the electronic properties of MnSi^ as can be seen in 
the manganese L2.3 x-ray absorption and XMCD spec- 
tra. These probe local 2p — > 3d transitions at the man- 
ganese dopant and therefore reflect its electronic struc- 
ture and magnetic moment. In Fig. [TJ exohedral clusters 
with n = 7 — 10 show nearly identical x-ray absorption 
spectra that indicate a very similar electronic structure of 
the manganese dopant. In contrast, the x-ray absorption 
spectrum and thus the local electronic structure is more 
complex and varies strongly with the number of silicon 
atoms for endohedrally doped MnSi^ . Yet more striking, 
exohedral MnSi^ shows a pronounced XMCD asymme- 
try that vanishes for endohedral species. Even without 
applying XMCD sum rules [2T] , the XMCD asymmetry is 
a qualitative and direct probe of magnetism and clearly 
indicates that manganese in exohedrally doped silicon 
clusters carries a magnetic moment which is quenched 
upon encapsulation. The residual XMCD asymmetry, 
which is observed for n = 11 and 12, is assigned to a 
slight contamination with < 5 % of Mn2Si„_ 2 , because 
this XMCD signal and the corresponding lines in the x- 
ray absorption spectrum were tested to be proportional 

that was observed simulta- 
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FIG. 2. (color online) Eigenvalue spectrum (bars), total 
(dotted line), and local Mn 3d-projected (solid line) DOS with 
isosurface plots (at ± 0.04 a.u.) of the Mn 3d orbitals and 
the highest occupied orbital for MnSi^~ and MnSi^. Positive 
(negative) values represent the spin-up (spin-down) channel. 
The d-like orbitals are localized on the manganese atom in 
MnSiy , whereas they are delocalized in MnSif 4 . 



To further analyze the exohedral-to-endohedral and 
magnetic-to-nonmagnetic transition in MnSi^, we per- 
formed a thorough global geometry optimization in a 
simulated annealing [53] and modified "big bang" [55] 
approach, details of which are given in the Supplemen- 
tal Material [26]. The calculations were carried out in a 
DFT framework using the Perdew-Burke-Ernzcrhof one- 
parameter hybrid (PBE0) [57] as implemented in tur- 
BOMOLE [5H]. The PBE0 exchange correlation (xc) func- 
tional was chosen because it partly cancels the effects 
of the self-interaction error [29] that is inherent in com- 
monly used semilocal functionals and often leads to er- 
roneous results for the electronic structure of systems in 
which the highest occupied orbitals differ significantly in 
their degree of spatial localization [3U]. This is partic- 
ularly true for systems containing transition metal ele- 
ments, such as MnSi+. 

The assigned ground state structures of MnSi^ that re- 
sult from our calculations are depicted in Fig. [T] With 
the exception of MnSi^~, the predicted structures of ex- 
ohedrally doped MnSi^ clusters for n = 7—10 corre- 
spond to those of Si+ +1 (cf. Refs.H9land[3T]). where one 
silicon atom is replaced by manganese. In contrast, no 
structural similarity to the corresponding Si^ +1 clusters 
can be observed in the endohedral size regime, where 
manganese is encapsulated by silicon. These structural 
findings agree qualitatively with results of the DFT and 

[13] . More- 



infrared spectroscopy study by Ngan et al. 



to the amount of M^Si^ 
neously in mass spectrometry. 



over, our calculations show that the magnetic moment of 
MnSi^ is quenched from 4 fis to /is at the exohedral- 
to-endohedral transition, in perfect agreement with the 
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FIG. 3. (color online) Eigenvalue spectrum (bars), total DOS 
(dotted line), and manganese-projected local DOS (solid line) 
of the true endohedral ground state of MnSi^ (top) and of 
the false exohedral ground state predicted by PBEO (bottom). 
Aei is a measure for how much the PBEO eigenvalues are 
affected by the self-interaction error (see text). 



FIG. 4. (color online) Magnetic moment versus weighted co- 
ordination doN c /a for ground state (open circles) and higher 
energy (solid circles) isomers of MnSiJ (n = 7 — 14), isolated 
neutral Mn impurities in bulk ([4, solid squares) and amor- 
phous silicon ([BJ, open triangles), and in a silicon nanocrystal 
([8], solid diamond). Inset: bond length distribution of MnSig" 
with first coordination sphere (dotted line). 



experimentally observed disappearance of the XMCD sig- 
nal from n < 10 to n > 11. 

The change in local electronic structure of the manganese 
dopant that was observed in the x-ray absorption spectra 
at the structural transition of MnSi^ in Figfllis reflected 
in the occupied eigenvalue spectrum, which is shown in 
its usual interpretation as a density of states (DOS) in 
Fig. [2] for MnSiy~ and MnSiJj, representing the exohedral 
and the endohedral size regime. For MnSi^, the occupied 
states of manganese 3c? character are mostly isolated at 
ps 4 — 5 eV below the Fermi level and are tightly local- 
ized at the manganese dopant, i.e., they largely preserve 
a (perturbed) atomic character and only weakly interact 
with silicon states, as can be seen from the DOS and the 
isosurface plots. Because of this weak interaction, the 
localized orbitals of manganese 3c? character are qualita- 
tively very similar for all exohedral MnSi^" clusters and 
lead to the nearly identical x-ray absorption spectra in 
Fig. [l] This notion of, at least partly, atomic manganese 
3c? states is lost in the endohedral size regime, repre- 
sented by MnSi^ in Fig. [2] Here, orbitals with partial 
manganese 3c? character are strongly hybridized with sil- 
icon states and are shifted to ps 0.5 — 2 eV below the 
Fermi level, i.e., they participate strongly in bonding 
and are delocalized over the silicon frame. Consequently, 
the manganese 3c? derived partial DOS sensitively de- 
pends on the structure of the silicon cage, which is re- 
flected in the variation of the x-ray absorption spectra 
of endohedral MnSi^ for different n in Fig. [I] As a re- 
sult of strong spd hybridization with the participation of 
all manganese valence orbitals, the magnetic moment is 
completely quenched in endohedrally doped clusters. 
Because of the well known self- interaction error [3§], a 
treatment of this change from a partly localized, atomic- 



like situation to one in which all orbitals are delocalized 
on similar length scales poses serious difficulties to stan- 
dard approximations of the xc functional. Therefore, the 
DFT results have to be analyzed carefully at the struc- 
tural transition around MnSi ^ , where exohedral and en- 
dohedral isomers can be expected to be closest in en- 
ergy. For MnSi^, PBEO predicts an exohedral ground 
state with a total magnetic moment of 2 /i^, which is 
1.13 eV lower in energy than an endohedral isomer that 
would be in agreement with the experimental results. We 
attribute this discrepancy to uncertainties in theory for 
two reasons: First, we argue that the experimentally ob- 
served endohedral structure is indeed the ground state 
since we have observed ground-state structures of the re- 
lated systems Si^ and VSi^ [TSl E3 H2] under similar 
experimental conditions and for comparable size ranges. 
Second, a closer look at the theoretical result reveals that 
the failure of PBEO for MnSi^ can indeed be explained 
in terms of the effect of self-interaction on the eigenvalue 
spectrum of both isomers. The self-interaction error a 
of the orbital ifi, 

e i = (c Pi \v H [|^| 2 ] +<r° x [l^| 2 ,0]|^), 

can be used to quantify the reliability of the eigen- 
value spectrum [30] . Here vh is the electrostatic Hartree- 
potential and v^ prox is an approximate xc potential. De- 
tails of the calculation of ej are given in the Supple- 
mental Material and references therein [33 -35 . The 
self-interaction corrected eigenvalues &j can be estimated 
as Si ps e ^ pprox — ei [29], where £^ pprox results from a 
self-consistent calculation using w^ pprox . Fig. [3] compares 
the difference Aei between Ei and the PBEO eigenvalues 
£ pbeo ag a measure f now reliably PBEO cancels the 
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orbital self-interaction error in the two MnSij^ isomers. 
The large value of Ae, as well as its scatter shows that 
the exohedral PBEO result is nonreliable in the case of 
MnSi^ because self-interaction strongly affects orbitals 
that participate in bonding. For the endohedral isomer 
with delocalized orbitals, the PBEO eigenvalues are al- 
most identical to £j, and are thus reliable. In view of 
these arguments, we conclude that MnSi^ is the small- 
est endohedral structure. This reassessment of the ener- 
getic order does not alter the observed quenching of the 
magnetic moment at the exohedral-to-endohedral tran- 
sition since every endohedral isomer of MnSi„ that re- 
sulted from our calculations is predicted to adopt a sin- 
glet state by PBEO. We therefore stress that these delib- 
erations only lead us to reconsider the predicted critical 
size in accordance with our reactivity studies, but do not 
change the electronic or magnetic properties. 
The structural change at the exohedral-to-endohedral 
transition of MnSi^ can be quantified by the coordina- 
tion number N c of manganese, i.e., the number of sil- 
icon atoms in the first coordination sphere as exempli- 
fied for MnSi^t" in the inset of Fig. [4] In exohedral clus- 
ters, manganese adopts a minimal coordination number 
of N c — 2—4, while in endohedral clusters N c is maxi- 
mized to 11 — 14, i.e., all silicon atoms are within the first 
coordination sphere of manganese. The average Mn-Si 
bond length a elucidates why encapsulation of the man- 
ganese dopant becomes energetically favorable only for 
n > 11: In the ground state structures, the Mn-Si near- 
est neighbor distance expands from a — 2.43 — 2.58 A in 
exohedral to a — 2.53 — 2.67 A in endohedral clusters. In 
contrast, it would be compressed to a = 2.35 A in the 
higher energy endohedral isomer of MnSi^~ . Even though 
manganese favors high coordination in silicon |36j . this 
strain, which becomes even more pronounced in smaller 
clusters, precludes endohedral ground states for n < 10. 
The abrupt change in coordination at the structural tran- 
sition is interrelated with the quenching of the mag- 
netic moment as illustrated in Fig. [4j Here, the calcu- 
lated magnetic moments of MnSi^ are plotted versus the 
weighted coordination number doN c /a, which takes into 
account both the number of silicon nearest neighbors iV c 
as well as their average distance a to the manganese 
atom, normalized to the nearest neighbor distance do 
in bulk silicon. Low-coordinated exohedral clusters with 
doN c /a — 1.9—3.7 (N c — 2—4) carry a magnetic moment 
of 4 /xb, which is quenched to /xb m high-coordinated 
species with d„N c /a = 10.2 - 12.3 (N c = 11 - 14). This 
relation of magnetic moment and weighted coordination 
also holds for higher-energy isomers that are included in 
Fig. [4] and mark the transition from magnetic to nonmag- 
netic impurities around doN c /a ~ 4. 
Fig. [4] shows that this observation can be generalized 
to extended systems, i.e., to a neutral manganese im- 
purity at a substitutional site in crystalline silicon ^ 
or in hydrogen-passivated silicon nanocrystals [8] . It also 



applies to very low concentrations of manganese in amor- 
phous silicon, for which a possible relation between mag- 
netic moment and coordination has been pointed out [6] . 
However, although N c is the leading term, it does not ac- 
count for the dependence of the local magnetic moment 
on the nearest-neighbor distance |3] that becomes impor- 
tant around N c — 4 and is included in doN c /a. As can be 
seen in Fig. [4j manganese-doped bulk silicon is just at the 
transition from high-spin to low-spin states and there- 
fore reacts very sensitively to changes in doN c /a. This 
might explain the large scatter in experimental results 
on manganese-doped silicon and indicates that high-spin 
states could be stabilized by an appropriate expansion 
of the lattice parameter, e.g., in ultrathin films or passi- 
vated nanocrystals. 

In summary, the magnetic moment of manganese-doped 
silicon has been investigated over a wide range of struc- 
tural parameters, including extreme coordination num- 
bers from 2 - 14. The study of singly doped, size-selected 
MnSi^ clusters avoids impurity-band formation or inter- 
action between impurities that might be present in ex- 
periments on bulk samples, but also in calculations with 
periodic boundary conditions. We are thus able to show 
that the observed quenching of the magnetic moment is 
not a result of impurity band formation but of the elec- 
tronic interaction with the silicon host. A universal cor- 
relation of the magnetic moment and the weighted co- 
ordination number is observed, providing guidelines to 
the stabilization of high-spin states in dilute manganese- 
doped silicon. 
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